Singularity formation and blowup of 
complex- valued solutions of the modified KdV 

equation 
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Abstract. The dynamics of the poles of the two-soliton solutions of the 
modified Korteweg-de Vries equation 

u t + 6u 2 u x + u xxx = 

are determined. A consequence of this study is the existence of classes of 
smooth, complex-valued solutions of this equation, defined for — oo < x < oo, 
exponentially decreasing to zero as \x\ — > oo, that blow up in finite time. 



1 Introduction 

Studied here is the modified Korteweg-de Vries equation 

u t + 6u 2 u x + u xxx = 0, (1.1) 

which has been derived as a rudimentary model for wave propagation in 
a number of different physical contexts. The present paper is a sequel to 
the recent work [9] wherein the dynamics of the complex singularities of the 
two-soliton solution of the Korteweg de Vries equation, 

u t + 6uu x + u xxx = 0, (1.2) 
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were examined in detail. 

The study of the pole dynamics of solutions of the Korteweg-de Vries 
equation and its near relatives began with some remarks of Kruskal [H] in 
the early 1970's. More comprehensive work was carried out later, see |18j . 
[TT] and [12]. One goal in our preceding paper [9] was to understand in more 
detail the propagation of solitons in a neighborhood of the interaction time. 
Another motivation was an idea to be explained presently concerning singu- 
larity formation in nonlinear, dispersive wave equations. More particularly, 
we are interested in both the generalized Korteweg-de Vries equation 

u t + u xxx + u p u x = (1.3) 

and coupled systems of Korteweg-de Vries type, viz. 

f u t + u xxx + P(u, v) x = 0, ^ , 

\ v t + v xxx + Q(u,v) x = 0, 

where P and Q are, say, homogeneous polynomials. 

Concerning singularity formation, it is an open question whether or not 
smooth, rapidly decaying, real-valued solutions of (11.31) develop singularities 
in finite time in the supercritical case p > 5. In the critical case p = 4, blowup 
in finite time has been established by Martel and Merle [17], [16] whilst for 
subcritical values p = 1,2, 3, there is no singularity formation for data that 
lies at least in the Sobolev space if^IR). (However, solutions corresponding to 
infinitely smooth initial values lying only in L 2 (1R) can develop singularities; 
see [7J.) Numerical simulations reported in [I] of solutions of (\1.3\i initiated 
with an amplitude-modified solitary wave indicate blowup in finite time. 
Such initial data has an analytic extension to a strip symmetric about the 
real axis in the complex plane. The results of Bona, Grujic and Kalisch [13] , 
[5] and [5] indicate that blowup at time t has to be accompanied by the width 
of the strip of analyticity shrinking to zero at the same time. This points 
to the prospect of a pair of complex conjugate singular points coliding at 
a spatial point on the real axis, thereby producing non-smooth behavior of 
the real-valued solution. It was shown in [9] that, in certain cases, curves of 
singularities do merge together. This happens at the moment of interaction 
of the two solitons when the amplitudes are related in a particular way. 
This result provides some indication that the blowup seen in the numerical 
simulations might occur because of the coalescence of curves of complex 
singularities. Such ruminations seem to justify interest in the pole dynamics 
in the context of (jl.3p . even in the case where the initial data is real-valued. 
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If one considers instead complex-valued solutions, it was shown in jH] 
(and see also [TU]) that, in the case of spatially periodic boundary conditions, 
equation (II .3p has solutions which blow up in finite time for all integers p > 1. 
Explicit examples of smooth, complex-valued solutions defined for x,t E M. 
of the Korteweg-de Vries equation (jl.2p that blow up in finite time have 
been given in 0, [El EH] • One outcome of the present paper is an explicit 
example of a blowing-up soution of the modified Korteweg-de Vries equation 

ALU). 

For the system (ll.4p where P and Q are homogeneous quadratic poly- 
nomials, conditions on the coefficients are known that imply global well- 
posedness for real-valued initial data (u ,Vq) (see [5]). And, the pole dy- 
namics investigated in [9] for the Korteweg-de Vries equation (II. 2p itself, 
(KdV-equation henceforth) reveals that the choice 

P(u,v) = u 2 — v 2 and Q(u,v) = 2uv (1.5) 

leads to a system (11.41) possessing solutions that develop singularities in finite 
time. (This is the system that obtains if complex-valued solutions of the 
KdV-equation are broken up into real and imaginary parts.) 

The latter result is obtained by a careful study of the pole dynamics of the 
explicit two-soliton solution U = U(z, t) = U(x + iy, t) of the KdV-equation 
in the complex z-plane. It transpires that as a function of time, most of 
the singularities of this exact solution, which are all poles, move vertically 
in the ^-direction in the complex plane as well as propagating horizontally 
in the ^-direction. As a consequence, by choosing y appropriately, one can 
arrange that the function u(x,t) = U(x + iyo,t) is a complex-valued solution 
of the KdV-equation that, at t = 0, is infinitely smooth and decays to zero 
exponentially rapidly as x — > ±oo, but which blows up for a positive value 
t > 0. If we write u(x,t) = v(x,t) + iw(x,t), then the pair (v,w) is a solution 
of (II .4p with the choice (I1.5P that starts at (t> , w ) = {v (•, 0), w(-, 0)) smooth 
and rapidly decaying, but which forms a singularity in finite time. It is worth 
noting that by an appropriate choice of the particular two-soliton solution, 
the initial data (vo,wq) can be taken to be arbitrarily small in any of the 
usual function spaces used in the analysis of such equations. 

The motivation for the current paper was to see if the phenomenon just 
described, i.e. curves of singularities merging together, can occur for the 
modified KdV-equation, and if this behavior might provide additional insight 
into the possible ways a singularity is produced in nonlinear, dispersve wave 



3 



equaions. Indeed, it turns out that the same phenomenon does occur. In 
certain cases, curves of singularities do converge together at one point. On 
the other hand, the behavior of these curves is not so different from what 
occurs in the context of the KdV-equation. 



2 The two— soliton solutions 

Preliminary analysis of the two-soliton solutions of the mKdV-equation 

Ut + u xxx + 6u 2 u x = 0, x G C, t E R, (2.1) 

are set forth here in preparation for the investigation of their pole dynamics. 
We begin with a standard transformation enabling one to express solutions 
of (12. ip in terms of solutions of another equation. Let u = v x where u is a 
solution of ( 12. ip . Then v satisfies the equation 

^ (vt + v xxx + 2vl) = 0. 

Assume v and its derivatives vanish at infinity and search for solutions of 
the latter equation of the form v = 2arctan(g). A calculation shows that v 
satisfies 

vt + v xxx + 2v 3 x = (2.2) 

if and only if 

(i + g 2 )(gt + 9xxx) + Qgx(g 2 x - gg X x) = o. (2.3) 

This yields a solution to (12. ip having the form 

u(x,t) = 2( a rct a ng(x,t)) x = (2.4) 

/x l + g(x,ty 

It is important to note that equation (12. ip and (12.31) are both invariant under 
change of sign. That is, u is a solution of (12. ip if and only if — u is a solution 
and, likewise, g is a solution of (I2.3P if and only if —g is a solution. Also, 
replacing g by 1/g in (12. 4 p has the effect of multiplying the solution u by —1. 
More precisely, if u is given by (12.41) . then 

2(arctan— ^— ) = -2 „ 9x ( X,t \ = - u (x,t). (2.5) 
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The well-known soliton solution of ( 12.1 j) has a hyperbolic secant profile. 
In detail, for any amplitude value k > 0, it is straightforward to check that 

g(x, t) = exp ( — k(x — x ) + k z t) = exp ( — k(x — x — k 2 t)) (2.6) 

is a solution of ( 12. 3ft . The corresponding solution u of ( 12. ip is the soliton 
solution with speed /c 2 and is given explicitly as 

g— k(x— xo)+k 3 t 

u(x,t) = 2(arctan^(x,t)) ;c = -2k - - e „ 2k{x „ Xo)+2kH (2.7) 
= -ksech(-k{x-x ) + k 3 t). (2.8) 

Note that the choice of sign in the exponential function g in (I2.6P produces 
the negative soliton. Replacing g by either — g ovl/g will produce the positive 
soliton. 

The two-soliton solutions are a little more complicated. The formulation 
presented here is based on that appearing in [I]. As just noted, there are 
both positive and negative soliton solutions of ( 12. ip . Consequently, there are 
two types of two-soliton solutions, namely interacting solitons of the same 
or opposite signs. 

Suppose first that < k\ < fc 2 . Define the functions fj by 

f){x, t) = exp{-k jX + kp), j = 1, 2. (2.9) 

Of course, this definition omits two arbitary spatial translations; these will 
be added later. Define two auxilliary functions, g + and g~ by 

where 

k 2 - k\ 

Proposition 2.1. The functions g + and g~ defined in Ii2.10\) - f2.11\) are so- 
lutions to A2.3\) . 
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Proof. It suffices to provide the proof for g~. Indeed, if one replaces f x by 
— fx, then g~ is transformed into g + , and all the calculations below remain 
valid in this case. 

Temporarily, set g = g~ . Notice that fj X = —kjfj and fjt = k^fj for 
j = 1,2. Thus, the quantities g t , g x , g X x and g xxx may be expressed in 
terms of fj and kj. First, differentiate with respect to time and come to the 
expression 

k\h ~ kjh + kjhfl - kjfih 
9t 7 (I + /1/2) 2 

Similarly, the derivative with respect to x is 

k 2 f2 - hh - hfJl + k 2 flf 2 
9x 7 (I + /1/2) 2 

Differentiating the latter expression leads to 

9" = n TTTYs te/i-^V 2 +(A: 2 +4fc 1 A: 2 +fc 2 2 )(/j|-/ 2 / 2 )-A: 2 A 2 /|+fc 2 2 / 1 3 / 2 2 
Differentiating once more and simplifying gives 
9*** = n -TTY4 N/ 2 - k\h - {k\ + 4k 3 2 + Qk\k 2 + I2k x k 2 2 ){f x f 2 + /j 3 / 2 

U + J1J2) L 

+ (4fc? + k\ + §k x k\ + I2k\k 2 ){f x f 2 + flfl) - k\flf A 2 + klftfl 
It follows that 

9t + 9m = V^t^ [hf 2 j2 ~ hhfl + hflf 2 - hfffl] 
U + J1J2) 

on the one hand, and 

9l - 99 XX = n 2) f u ~ k *) 2 (fih + f!fi)- 

U + J1J2) 

Shhfif* + (h + k 2 )\flf 2 + hfl 

on the other. It is now straightforward to check that g satisfies equation 

□ 
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The above proposition implies that 



m ± (x, t) = 2( arctan g (x, t)) x = — 




are solutions to (12. ip . A further calculations reveals that 




(2.13) 



and 




(2.14) 



where the new combinations 



G + = A; 1 / 1 (l + / 2 2 ) + A ;2 / 2 (l + / 1 2 ), 
G- = -k 1 f 1 (l + f^) + k 2 f 2 (l + f!), 



(2.15) 
(2.16) 
(2.17) 
(2.18) 



^ + = (l-/i/ 2 ) 2 + 7 2 (/i + / 2 ) 2 , 

F- = (i + A/ 2 ) 2 + 7 2 (/i-/ 2 ) 2 , 



have been introduced. Note that the functions with a superscript "+" are 
obtained from the functions with a superscript "— " simpy by replacing fx 
by —fx- If every occurence of kx is replaced by —kx in formula ( 12.141) for 
u~, then fx is replaced by 1/fx and 7 is replaced by I/7. Simplifying the 
resulting expression gives exactly the formula (12 . 1 3[) for u + . In other words, 
the formulas for u + and u~ can be obtained from each other by replacing 
every occurence of kx by —kx- 

The function u + is a two-soliton solution of (12.1 1) having two positive 
interacting solitons whist u~ is a two-soliton solution with two interacting 
solitons of opposite sign, the faster one being the positive one. One can see 
this graphically using MAPLE or Mathematica, for example. Analytically, 
there is an explicit relationship between the fomulas for given by (12.131) 
and ( 12.1411 and the single soliton solutions of speeds k\ and k\. As in (12.7j) . 
let Uj be the centered, positive, soliton solution 



of speed kj,j = 1,2. If the formulas for (12. 13ft and (12. 1411 are both divided 
by (1 + fl)(l + /|), there obtains 



Uj(x, t) = kj sech(—kjX + kj 



i + fj(x,ty 



(2.19) 



Ml + u 2 



(2.20) 



u =7 



D+ 



7 



and 

-u 1 + u 2 

u =7 — — (2.21) 

where 

(l-/i/ 2 ) 2 + 7 2 (/i + / 2 ) 2 



D~ 
D~ 



(l + /i 2 )(l + /|) 
(l + / 1 / a ) 2 + 7 3 (/i-/a) s 



/i 2 )(i+/i: 



Formulas (I2.20p and ( I2.2ip show in particular that u + (x, t) > for all x G K 
and that m~(x, t) > precisely for those a; G K and £ G R for which ^(x, £) > 
ui(x, £). 

In Section HI the asymptotic behavior of the singularities of are exam- 
ined for large positive and negative time. We will see that they separate into 
two groups, corresponding to the two single solitons. More remarks on the 
shape of the two-soliton solutions during their interaction are to be found in 
Section [SJ 

Before ending this section, we return to the issue of spatial shifts. For 
< ki < k 2 and X\, x 2 G R, let 

fi(x,t) = exp(—ki(x - aci) + k\t), 
f 2 (x, t) = exp(-k 2 (x - x 2 ) + k^t), 

and 

h(x,t) + f 2 {x,t) 



g+(x,t) = 
g~(x,t) 



l- 7 - 2 /i(M)/2(x,t) 
fi(x,t) - f 2 (x,t) 



l + ^- 2 fi(x,t)f 2 (x,t) 
Define the interaction time to and the interaction center xq of g ± (x,t) to be 

io = -|3~feTIoW og7 a,ld (2 ' 22) 
x » = 14-ki ~ (k 2 + w5 log7 ' (2 ' 23) 

respectively. 
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Proposition 2.2^Let t andxo be the interaction time and interaction center 
for g + (x,t) and g~(x,t). Then, for any t G M. and x G C, we have 
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(x,t) = g ± (x -x ,t- t ), 



(2.24) 



where g^ 1 are defined in Ii2.10\) - f2.11\) . Moreover, the functions M ± (-,to) — 
2(arctang ± (-, t )) x are symmetric about the point xq on both R and C. 

Proof. It suffices to find (xo, t ) G M 2 such that fj(x, t) = jfj(x — xo, t — t ) 
(j = 1, 2) for all t G R, x G C. Equivalently, fl2^4"|) will be satisfied if 



Since (xo,t ) given by fl2.22p - fl2.23p is the solution for this system, the first 
assertion is proved. To see the symmetry of the functions M ± (-,t ) about x , 
it is only necessary to deduce from (123311 that m ± (-,0) = 2(arctan(^ ± (-, 0))* 
is an even function. This is easily verified since 5 f± (— x,0) = —g ± (x,0) for 



Since (mKdV) is invariant under time- and space-translation, it is con- 
cluded from Proposition 12. 21 that the functions u ± (x,t) = 2(arctan(g ± (a;, t)) x 
are also solutions to (12. ip . The point is that, for a general two-soliton solu- 
tion, the time and place of the interaction are given by t and Xq in (12.221) 
and (12.231) . respectively. In particular, the solutions are already normal- 
ized so that the interaction time is t — 0, and the center of the interaction 
is at x = 0. It is interesting to note that the values of to an d Xq in the 
above proposition are precisely the same as for the two-soliton solution of 
the Korteweg deVries equation, as given in Theorem 1 of [9]. 

3 Singularities of the two-soliton solutions in 
the complex plane 

As mentioned in the introduction, the two-soliton solutions u(x,t) of (12. ip 
are viewed as meromorphic functions in the complex variable x. We aim to 
determine how the dynamics of the singularities in C reflect the behavior of 
u(x, t) when x is restricted to the real axis. 




klXo—kf 
k2XQ—k2 




all x G C. 



□ 
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Consider first the (one)-soliton solution, given by (I2.7p . It is immediate 
that the singularities of these solutions are precisely 



x = x + k 2 t + — - (3.1) 

where m runs through the odd integers. Moreover, these singularities of u 
are all simple poles. Thus, the speed of the soliton is exactly the speed of its 
poles in the complex plane, while the position of the maximum point of the 
soliton at time t is the real part of the position of these poles. The imaginary 
part of the singularity remains constant in time of course. 

The poles of the two-solitons solutions ( 12 . 13f) and ( I2.14p . correspond to 
the zeros of the functions F ± and G ± defined in fl2~T5|) . fl2~T6|) . ( TXTTj) and 
(I2.18p . It will turn out that the singularities of u^ 1 are all simple poles, just 
as for the one-solitons. In all but a specfic class of exceptional cases, these 
poles correspond to simple zeros of F^. 

Recall that ki and k 2 are called commensurable if there exist positive 
integers pi and p 2 such that 

T- = -■ ( 3 - 2 ) 

Without loss of generality, we may take it that gcd(pi,£>2) = 1- In what 
follows, we will always assume that the integers p\ and p 2 appearing in (13. 2p 
are without common prime factors. In this case, F (x,t) and G (x, t) (and 
thus M ± (x,t)) are periodic in x with minimal imaginary period 2n\i where 

* = M- < 33 > 

Lemma 3.1. For any tGl, the zeros of F ± (-,t) are simple, except for the 
following special case. If k\ and k 2 are commensurable, and if Pi,p% G N and 
A > are given by $3. S\) and {i3.3\) . with pi ± p 2 G 4M ; then there is a fourth 
order zero of F ± (-, 0) at x = (| + q)\ni for all q G Z. 

Proof. From (I2.18p . we infer 

F- = 2[(1 + fif 2 )(fij2 + hf2 X ) + 7 2 (/i - h)U*> ~ /2x)]- 

Noticing that fj X = —kjfj, (j = 1,2) and setting X = fi and Y — f 2 , it 
follows that if x is a zero of F~(-, t) of order greater than or equal to 2, then 

f (1 + XF) 2 + 7 2 (X-F) 2 = and 

1 -(/=! + fc 2 )(l + Xy)XF + 7 2 (X - F)(A; 2 F - ^X) = 0. [ ° ] 
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Let (X, Y) be a solution of this system. Then, from the first equation we 
deduce that 

1 + XY = iei(X-Y), (3.5) 

where e is either 1 or -1. Inserting this into the second equation, it is found 
that 

ie(k 2 - k x )XY = k 2 Y - hX, 

provided X ^ Y. Extract the product XY and inject it in (13. 5p to obtain 
the linear relation 

Y = —X + IE : . 

k\ k\ 

Then formla H 3 . 5 1) implies that X must satisfy 

(k 2 v t . k 2 - kA . ( ( k 2 \ . k 2 -k x 
1 + XI - — X + %e — : = tej 1 — — } X — te 



k\ ki J \\kij ki 

which simplifies to 

X 2 + 2ieX -1 = 0. 
It follows that X = Y = ±i. Thus, it must be that e - k ^ x + k \ l = e -fca>+fc£t 
■ /'. from which we deduce that — k\x + k\t and — k 2 x + k 2 t are both purely 
imaginary. Therefore, t = and x is purely imaginary. Additionally, we have 

k\x = [\ + qi)iri, f „ 
, ){ , ; . for some qi, q 2 G L. 
k 2 x = (| + g2)7r?, 

Thus k\ and /c2 are commensurable and in fact, k 2 jk\ = r 2 jr\ with rj = 
1 + 2qj. Noticing that q 2 — q\ G 2Z (because fcia; — k 2 x G 2inZ), it transpires 
that r2 — ri = 4fc for some G Z. It follows that d = gcd(ri,r2) divides 
4k and therefore <i divides k since <i is necessarily an odd integer. Define 
Pi and p 2 by A^Ai = P2/P1 with gcd(pi,p 2 ) = 1- With this definition, 
P2 — Pi — {r 2 — ri)/d = Ak/d G 4N. There are precisely two such complex 
numbers in the fundamental strip S = {x G C, — Xtt < Imx < Xtc}, and they 
are x = ±^ni. 

It is straightforward to ascertain that for such values of x, F~ x (x,0) = 
F~ xx (x, 0) = and F~ xxx (x, 0) ^ 0. The same arguments hold for F + . □ 

Lemma 3.2. For any t G M, the zeros of F~(-,t) and G~(-,t) are distinct, 
except the following special case. If k\ and k 2 are commensurable, and if 
Pi,p 2 G N and A > are given by k3. || ) and k3. 3\) . with p\ ± p 2 G 4N, then 
there is a third order zero of F~(-,0) at x = (| + q)\ni for all q G Z. The 
same statement is true with F~ and G~ replaced by F + and G + . 
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Proof. In the notation of the proof of Lemma 13.11 the result follows if the 

system 

(i + xy) 2 + 7 2 (x-r) 2 = o, 
-hxn + Y 2 ) + k 2 Y(i + x 2 ) = o, {6 - 0> 



admits as its only solution X = Y = ±i. As before, if (X, Y) is a solution, 
there exists e G { — 1,1} for which 

1 + XY = ie'y(X-Y). 

Extracting from this the variable Y and computing 1 + Y 2 leads to 

Y - l£lX ~ 1 
X + ie-y 

and so 

l + y 2 = (l- 7 ^) 



i + x 2 



(X + ^ 7 ) 2 

Now insert this into the second equation of the system (13. 6p and, assuming 
by contradiction that 1 + X 2 is not zero, simplify the outcome. It follows 
that 

- ! 2 )X = (te^X - 1)(X + ie 7 ). 
After further simplifications, this becomes 

X 2 + 2ieX -1 = 0, 
and the claim follows. □ 

Definition 3.1. The values ofk\ and k 2 for which -F ± (-, 0) has multiple zeros 
are collectively referred to as the exceptional case. This occurs when k\ and 
k 2 are commensurable, pi and p 2 are odd, and 

p 2 - pi E 4N 

when considering F~ , and 

P2 + Pi e 4N 

when considering F + . 
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Lemmas 13. II and 13.21 show that all singularities of u (•, t) are simple poles. 
In all but the exceptional case, they correspond to simple zeros of F . In 
the exceptional case, they correspond to a third order zero of ^ ± coinciding 
with a fourth order zero of F^. It follows from Rouche's Theorem that in 
the exceptional case, four simple poles converge, as t — y 0, to the simple pole 
at (| + q)Xiri for all q G Z. Moreover, by the residue theorem, the sum of 
the residues of u ± (-,t) at these four poles converge, as t — > 0, to the residue 
of the simple pole at (~ + q)Xni for all q e Z. 

Proposition 3.1. In all but the exceptional case, the poles o/w ± (-,t), i.e. 
the zeros of F ± (-,t), are described by analytic curves x : K — > C. In the 
exceptional case, these curves are defined and analytic separately for t < 
and t > 0. 

Proof. Consider the case of u + and F + . Since all the zeros of F + (-,t) are 
simple, the implicit-function theorem shows that, for a fixed time to, a zero 
of F + (-,t ) can be locally and uniquely continued as an analytic curve x{t) 
such that F + (x(t),t) = 0. Such a curve x(t) can be continued as long as 
it remains in a bounded region of C. Thus, we need to show that \x(t)\ 
must stay bounded as long as t remains in a bounded interval of 1R. First, 
it follows from Proposition 16.31 in Section [6] below that the imaginary part of 
x(t) must remain bounded. Furthermore, if Rex(t) — > oo in finite time, then 
the equation F + (x(t),t) = implies 1 = 0. If Rex(t) — > —oo in finite time, 
then the equation ^f^^^fj^^ f )2 = also implies that 1 = 0. 
A similar argument works for u~ and F~ . 

□ 

This section is closed with a specific example of an element in the excep- 
tional case, namely k\ = 1, k 2 = 5 (calculations done with MAPLE). In this 
case, A = 1 and F~ and G~ may be rewritten as 

G-(x, t) = -e 251 V 1 + 5e 12 V + 5e 12 V - e l y and 

F-(x,t) = e 25 V 2 + le 250t y w - ^e 126t y 6 + ~e 2t y 2 + 1, 

where y = e~ x . Note that at time t = 0, F~ and G~ are symmetric polyno- 
mials in y, and that i and —i are third order zeros of G~ , and fourth order 
zeros of F~. All the other zeros are simple. We also have the decomposition 

. , G-(x,0) -2 2 y-4i y + 4i 

0) = 3 „_; ' = : : + + 



F (x, 0) y — i y + i 2y 2 — iy — 2 2y 2 + iy — 2 
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4 Large— time asymptotic behavior of the sin- 
gularities 

In this section we show that the poles of u^, that is the zeros of F ± , separate 
out into two groups, as t — > ±00, behaving asymptotically as poles of single 
solitons, of speeds k\ and respectively This, of course, reflects the fact 
that are in fact "two-soliton" solutions of (12. ip . 

The situation is nearly identical to that obtaining for the two-soliton 
solution for the Korteweg-de Vries equation (compare the following theorem 
with Theorem 2 in [H]). In particular, the backward shift of the slower wave, 
well-known in the case of the KdV-solitons, is also present for the modified 
KdV-equation. 

Theorem 4.1. The asymptotic behaviors as t — )■ ±00 of the curves x(t) 
of zeroes of F ± (x,t) whose existence was determined in Proposition VS. 1\ are 
completely described as follows. 

1. For every odd integer m G Z, there exists a unique curve x ms -(t) of 
zeros of F ± (-,t) such that 

! TXVK% 

x m , s -(t) = k^t+— log T + — - + o(l), 

as t —00. 

2. For every odd integer m G Z, there exists a unique curve x m ^ s + (t) of 
zeros of F ± (-,t) such that 

= ^i^ - T lo S7 - 171— + o(l), 
k\ zki 

as t —7- 00. 

3. For every odd integer n G Z ; there exists a unique curve x n j- (t) of 
zeros o/F ± (-,t) such that 

j TVK% 

as t —7- —00. 
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4- For every odd integer n G Z ; £/iere exists a unique curve x n j+(t) of 
zeros o/-F ± (-,t) such that 

I TlTCl 

Xnj+(t) = kit + — log7 - — + o(l), 

as t — » oo. 

To prove Theorem 14.11 it is necessary to study the zeros of F ± (-,t) with 
respect to frames of reference which move at the speed of each constituant 
soliton. As in Section 4 of [9], we set 

z = x - kft, (4.1) 

w = x- k^t, (4.2) 

r = exp(A; 2 (A| - k\)t), (4.3) 

s = ey^(k 1 (k 2 2 -k 2 l )t). (4.4) 



If 

and 
then 



H ± (z, r) = (1 T re' {kl+k2)z ) 2 + 7 2 (e" fcl2 ± re~ k2Z f (4.5) 
J±( Wj s ) = ( s T e ~(k 1+ k 2 ) w y + 7 2( e -fc lW ± se -^)2 ; ( 4 g) 



F^x, t) = H^z, r) = s^^iw, s). 

Solutions z(r) of H ± (-,r) which remain bounded in C as r — > and as 
r — )• oo correspond to zeros x(t) of F ± (-,t) asymptotically traveling at speed 
k\ as t — > ±oo. Likewise, solutions w(s) of / ± (-, s) which remain bounded in 
C as s — > and as s — >■ oo correspond to zeros x{t) of F ± (-, t) asymptotically 
traveling at speed k\ as t — > ±oo. The following is true of the curves z(r). 

Proposition 4.1. For every odd integer m G Z, taere exists a smooth curve 
z^(r) defined in some interval ofr>0, such that H ± (z^ l (r),r) = and 

i . . 1 , miri , . m-i 4fc 2 h iu , n / 

z±(r) = -log 7 + — T (-l) — ^|Z^2^ 7 e 2fcl r + o(r) (4.7) 

as r — > + . For every ocM integer n G Z ; taere exists a smooth curve w^(s) 
defined in some interval of s > 0, such that J ± (w^(s), s) = and 



= -y- log 7 + 777- t -1 2 72 — rr~ 

fc 2 2fc 2 k% — k( 



<(s) = -- log 7 + — =F (-l)^7^72 7- fel/fe2 e^ nm 3 + o(s) (4.8) 
as s -)■ + 
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Proof. The relation H ± (z,0) = is satisfied if and only if there exists an 
odd integer m G Z such that 

1 mni 

and similarly, / ± (w,0) = is equivalent to 

1 nni 
w = — — log 7 + — — 

k 2 2k 2 

for some odd integer n G Z. Applying the implicit function theorem, there 
exist smooth curves z^(r) and w^(s) defined in a neighborhood of z = 
1/ki \og r y + mni/2ki and wo = — l/k 2 logj + nni/2k 2 , respectively, such that 
H ± (z^(r),r) = 0, z^(0) = z and I(w^(s), s) = 0, w^(0) = wq. It remains to 
calculate {z^'ir) and (w^)'(r). Differentiating the equation H ± (z^ l (r),r) = 
with respect to r yields 

[-(A; 1 +A; 2 )(lTre- (fel+fc2)2 )re- (fcl+fe2)2 +7 2 (e- fclZ ±re- fe22 )(A; 2 re- fc22 ±A; 1 e- fcl2 )]/ 
where z = z^(r). Taking r = gives 



±A;i7 2 e -2 fcl 2 m (0) (2 ±y (0) = (7 2 _ 1)e -(/c 1+ fc 2 )2±(0) ? 



and using that 



e -0i+fc 2 )2±(o) = I(_]_) n ^ 1 7 -fe2Ai e -'!rr"'"' 



7 

leads to 

(^)'(0) = T(-l)^p%7- fc2/fcl e-^' 



The derivative (w^)'(O) is similarly calculated. Differentiating the equa- 
tion J ± (w^(s), s) = with respect to s leads to 

[±(A:i+£:2)(sTe~ (fel+fc2)w )e-^ 

= -(sf e -(*i+*2)«') T 7 2 ( e - felW ± se^^e - * 2 '". 

At s = 0, there obtains 

[-(^ + fc 2 ) e -2(*i+fc 2 )^(o) - fc l7 2 e - 2fcl,u » (0) ]«)'(0) = T(7 2 - l) e -( fe i +fe2 )^(°) 
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and since 

e -(ki+ki)w£(p) 

we conclude 

(«£)'(«) = T ( 



A proof of Theorem 14.11 is now readily available. The last proposition 
immediately gives all the curves of zeros of F ± (-,t) claimed in Theorem 14.11 
The only remaining issue is to prove that there are no additional zeros of 
F ± (-,t). If ki and k 2 are commensurable, this is straightforward. For each 
t 6 1, F ± (-,t) is a polynomial in e~ xX of degree 2(pi +P2), and thus must 
have precisely 2(pi + p 2 ) zeros in the complex plane. The zeros of F ± (-,t) 
described in Theorem 14.11 account for all of them, for large \t\, and so no 
other zeros can exist. 

If ki and k 2 are not commensurable, the desired result can be obtained 
by approximating ki and k 2 by sequences {ki}^ =1 and {k 2 }^ =1) which are 
commensurable . 

5 The nature of the singularity in the excep- 
tional case 

In this section a more detailed analysis is undertaken of the singularity of 
in the exceptional case (see Definition 13. 11) . As described just after this 
definition, the singularity at (| + q)\iii, for any q G Z, corresponds to a 
fourth order zero of F ± (-,0) and is approached by four simple zeros of 
as t — > 0. The goal is to understand the behavior of a smooth curve of zeros 
of F(-,t) in a neighborhood of such a fourth order zero. 

We claim it suffices to analyse the singularity at ^m, i.e. the case q = 0. 
Indeed, since F ± (-,t) is 2n\i periodic, it is enough to consider q = —1,0. 
Furthermore, F (x,t) = F ± (x,t), and so F ± (x(t),t) = if and only if 
F ± (x(t),t) = 0. Hence, only the case q = need be examined. 
Remark that 

F ± (x,t) = e - 2 ^ +k ^ x+2 ^ +k ^ t F ± (-x,-t). 



□ 
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From this, it is deduced that F ± (x(t), t) = if and only if F ± (— —t) = 0. 
In other words, if x{t) is a curve of zeros approaching |7rz as t /* 0, then 
—x{t) is a curve of zeros approaching ^iri as t \ 0. 

Theorem 5.1. Suppose we are in the exceptional case wherein k\ and k 2 are 
commensurable and the odd integers Pi,p2 G N and A > are as in A3. || ) and 
j[3. 3]) . In the case of F + , it is assumed that pi + P2 G 4N and in the case of 
F~ , it is presumed that p\ — p 2 G 4N. Let x{t) be a smooth curve, defined 
fort close to 0, such that F ± (x(t),t) = and x(t) ^ ^iri but x(t) — > ^ni as 
t — > 0. Then, either 

(x(t) - ^nif , s 

lim^ 2 -— = -Yl 1 5.1 

or 

x(t) — ^iri ,9,9 , 

lim — 2 — = k 2 ,+ kl 5.2 

t^o t 

Remark 5.1. The similarity of the result in Theorem \5.1\ with the "excep- 
tional case" for the two-soliton solution of sthe KdV-equation is quite strik- 
ing. The behavior of the curves described in I5.il is exactly the same as given 
by Proposition 4-9 in JMj. However, the exceptional case for the KdV-equation 
corresponds to pi being an odd integer, and p 2 being an even integer. For the 
mKdV -equation, the "exceptional pole" is located at \^i, rather than at Xni 
as it is for the KdV-equation. 

Also, how does one explain the behavior described by li5.S\) ? In the case 
of the two-soliton solution of the KdV-equation, if pi is even and p 2 is odd, 
there is a horizontally moving pole approaching the singularity at n\i exactly 
as described by ( 15. This calculation was not carried out in J^j, but can 
be obtained by a simple modification of the proof of Proposition 4-9 in 
Thus, it appears that the exceptional case for the two-soliton solution of the 
mKdV-equation \2. 1\) includes the behavior of curves of singularities from 
two different cases of the two-soliton solution of the KdV-equation, namely 
the "exceptional case" , where p\ is odd and p 2 is even, as well as the case 
where p\ is even and p 2 is odd. It is precisely in these two cases that there is 
a pole located at the same place n\i (at time 0). 

Proof. The proof is provided for F~ . Similar arguments apply for F + . Let 
z(r) = x(t) — kft where r is defined in (14. 3p . It follows that z(r) is a smooth 
curve such that H~(z(r), r) = 0, where H is as in f)4.5p . Notice that z(r) — > 
^ni as r — > 1. Now, 

H~(z, r) = & 1 + re~ {kl+k2)z = ±2 7 (e~ fel " - re~ k2Z ). (5.3) 
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If the minus sign is chosen on the right-hand side of ( 15. 3p . we come to 

l^ e -kiz _|_ i 
T = . 

e -(fci+fe2)^ i^/e~ k2Z ' 

Differentiating this relation with respect to r, it follows that 

1 = [i^e^e-^ - rye"**) 

+ {i ie - klZ + l)(-(h + k 2 )e- {kl+k2)z + i^e-^)] z'(r). 



This may be rewritten as 

2^1 



i~fk 2 e- k2Z (1 + ie~ klZ ) . 
1 = —z , 



or what is the same, 



^g fc 22 ( e -(ki+k 2 )z _ - e -k 2 z\ 2 

z = i — '—. (5.4) 

fc 2 7(l + ie~ klZ ) 

Since pi and p 2 are both odd with p 2 —p\ £ 4N and z{r) — > ^ni as r — > 1, 
it must be that e~ klZ — >■ e~ P17 ™/ 2 and e~ fc22 — > e~ P27 ™/ 2 as r — V 1. If pi and p 2 
are both in 4N+ 1, then e _fcl2 and e _fcl2 both converge to — i as r — >■ 1, while 
if pi and P2 are both in 4N + 3, then e~ klZ and e~ fclZ both converge to i as 
r — )■ 1. We suppose first that p\ and P2 are both in 4N + 1. In this case, it 
follows from H5.4I) that 



( 7 + I) 2 
4/c 2 7 fc I _ k 't 



Ar) -> = 7^2 (5-5) 



as r — >■ 1, and thus 



z(r) - fvrz A; 2 

hm 



r-+i r - 1 k\-k[ 

Turning back to x(t), and using the fact that (r — l)/t — )■ ^2(^2 — ^1) as 
t — 7- 0, it follows that 

x(t) - £; 2 t - ^ni , „ 
hm — - — = ki 

t^o t 2 
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from which it is concluded that 



Um x(t)-| ! r ! = 

t->0 t 12 

If instead the plus sign is chosen on the right-hand side of (15. 3p . then it 
is immediately inferred that 



-k\z 



1 



e ~{ki+k 2 )z _|_ ^ e ~k 2 z ' 

Differentiating this equation with respect to r, we get 

1 

(^ e -{ki+k 2 )z _|_ ^ e -k 2 z^2 

-(i^e- klZ - 1)( - {h + k 2 )e- (kl+k2)z - i-fk 2 e- k2Z ) 



1 



which may be simplified to 

1 = - 



i~/k 2 e- k2Z (l -ie- klZ ) 2 . 

y 

(^ e -{k\+k 2 )z _|_ iry e -k 2 z^2 



From this, it is inferred that 



;i -%e-^ z fz' 



« '" 1Z ( e -(k 1+k2 )z + i7e -* a *)2 



hi 



Since pi and p 2 are both in 4N + 1, it is concluded that 



1 — ze 



2 — > 



k-2n k 2 {k\ - fc?) 



as r — 7- 1. A consequence of this is that 
d 



1 - ze 



-k\z\3 



3(1 - ze^ 12 ) 2 ^/^-* 12 )^ -> 



kiz\ „/ 



-12fc? 



L'Hopital's rule comes to the rescue and it is found that 

(l_ ie -fci*)3 _ 12 ^3 



lim 



>i r — 1 



- kl 



(5.6) 



(5.7) 
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Since 



lim ~~t — —i lim ^ = — the kl ^ m = k\, 



•-7/1/ - ^ — f^7ri Z 2^"^ 



we must have 



(z(r) - |7ri) 3 -12 

lim 



r->i r - 1 h{k\ - k\) ' 

Reverting to the original variable x(t) and using again that (r — l)/t — > 
h{k\ — k\) as t — y 0, it follows that 

t^o t 

and thus 

(x(t) - ^vu) 3 
lim — 2 —^- = -12. 

It remains to treat the situation where pi,P2 £ 4N + 3, in which case 
e -kjz(r) _s, ^ ag r _j, j n f ac ^ this case is "dual" to the one just treated and 
the same calculations lead to the result. To see this, assume first that the 
minus sign obtains on the right-hand side of (15. 3p . Then ( 15. 4p implies that 

(1 + «T fcl2 )V = -!^!(e-(*i+*»)* - %1 e~ k ^f -> (7 " 1)2 - " 4A;2 



hi hi h(kl - kf) 

as r — > 1. On the other hand, it is straightforward that 

1 + ie" fclZ , 
hm ^ — — = k\. 

z^^ni z — -^iri 

Following the same line of development as pursued for the positive sign in 
the previous situation where both pi and p 2 lie in 4N + 1 leads to 

(x(t) - ^tu) 3 
lim — 2 —^- = -12. 

t^-o t 

Now assume we have the positive sign on the right-hand side of (15.31) . From 
(15.61) we deduce 

, _ ie k2Z ( e - {kl+k2)z + i^e~ k2Z f _^ (7 + l) 2 _ h 



hi (l-ie- klZ ) 2 4hl k\-k\ 
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as r — > 1. Then, it is clear that 



lim 



x(t) 



/C-j^ ~\~ . 



□ 



6 Vertical movement of poles 

In this section we study the vertical motion of the singularities of the two- 
soliton solutions of ( 12. ip . As seen in Section [2], this comes down to 
studying the zeroes of F ± given by (I2.17p - (l2.18p . Recall the solution u + 
represents two interacting solitons of the same sign, while u~ represents two 
interacting solitons of opposite sign. Observe that 



and fi and f 2 are as in (12. 9p . Note the similarity in form between the above 
functions and the function F defined by formula (2.13) in [9]. 

Since F^~(x,t) = if and only if F 2 + (x,t) = 0, to study the zeros of F + , 
it suffices to study the zeros of F±. Similarly, since F{(x,i) = if and only 
if F 2 ~(x, t) = 0, to study the zeros of F~ , it suffices to study the zeros of Ff. 

As before, it is sometimes necessary to distinguish between the cases 
where k% and hi are commensurable and the cases where they are not. In 
the former case, denote by p\ and P2 the relatively prime positive integers 
such that (13. 2p holds, and let A be defined as in (13. 3p . Hence, if ki and k 2 
are commensurable, the functions f% and f 2 are 2^Xi periodic in x. Thus, in 
the commensurable case, all four of the functions F^, F 2 ± , must also be 27rAi 
periodic in x. 



F + (x,t) = F+(x,t)F 2 + (x,t) 
F~(x,t) = F^(x,t)F 2 ~(x,t) 



(6.1) 
(6.2) 



where 




1 + 17/1 + ^7/2 - /1/2 
1 - 17/1 - hh ~ /1/2 
1 + 17/1 - ^7/2 + fifc 
1 - «7/i + Hh + /1/2 



(6.3) 
(6.4) 
(6.5) 
(6.6) 
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Proposition 6.1. The zeroes of F + and F~ in the complex plane lie off the 
real axis. Moreover, in the commensurable case, if either F + (x,t) = or 
F~(x,t) = 0, then Imx ^ 2mii\ for all m G Z. 

Proof. Suppose x G E and F + (x, t) = 0. It follows that Ff(x, t) = F 2 + (x, t) = 
0. Taking real and imaginary parts produces the coupled system 

1 = fifi, 

fl — —f2, 

from which it follows that fi(x,t) 2 = —1, which is impossible. If x G M and 
F~~(x,t) = 0, it similarly follows that 

1 = -/1/2, 

fl — J2, 

from which one observes that fi(x,t) 2 = —1, which is again impossible. 

The last statement follows from the 27i\i periodicity which holds in the 
commensurable case. □ 

If k\ and &2 are commensurable, further information about the location 
of the zeros of F^ 1 can be obtained. 

Proposition 6.2. Suppose that k\ and k% are commensurable. If either 
F + (x, t) = or F~(x, t) = 0, then Imx 7^ (2m + l)7rA for any m G Z. 

Proof. Because of 27rAi-periodicity, it suffices to consider m = 0. Suppose 
Imx = 7rA and F^(x,t) = 0. Since x — x = 2rr\i, it follows from 2n\i- 
periodicity that F±(x,t) = 0. But this says precisely that F 2 + (x,t) = 0. 
Adding and subtracting the two equations, F^(x,t) = and F^~(x,t) = 
gives 

1 = /1/2, 

fl = ~f2i 

from which it is concluded that fi(x,t) 2 = —1, i.e. fi(x,t) = ±i and 
/ 2 (x,t) = =p. The latter system together with the definition ( 12. 9ft of the 
f/s implies 

fi,G oc — — k-^tj 
Pi,G X — — hj^t^ 
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and so t = Rex = 0. Thus x = n\i. However, 



/i(7rAz,0) = exp(— ki-n\i) = exp(— p^i) = ±1 



since p\ is an integer. This contradiction shows that F^~(x, t) can not be zero 
if Imx = 7rA. A similar argument applies to the other functions _F 2 + , F{ and 
F 2 . □ 

To show that most of the poles of the two-soliton solution of (12. ip move 
vertically, we need to reproduce calculations which are similar to those found 
in Section 3 of [9] . Unfortunately, it seems that the calculations in [9] do not 
directly apply to the present situation in all cases. In an effort at economy, 
we only prove here that, with certain very specific exceptions, the poles in 
the two-soliton solution always feature vertical movement. Consequently, 
we do not need to reproduce the entirety of Section 3 of [9J. Nonetheless, 
we cannot avoid certain, somewhat tedious calculations, closely modeled on 
those in [9j. 

The following notation 

a = — Imx, 

A kiRe x+k?t 

Til — e , 

— e , 

is taken from [0J. The function F± may be rewritten in this notation, viz. 

F+(x, t) = l + i-iA x e ikia + i^A 2 e ik2a - A 1 A 2 e i< - kl+k ^ a . (6.7) 

Attention is first focussed on the solution u + . To investigate possible ver- 
tical movement of poles of u + , we examine more closely the zeros of F^~(x, t). 

Proposition 6.3. Suppose F^(x,t) = 0. Then, coskia = if and only if 
cos k 2 a = 0. Moreover, the relation 

(^A 2 — — ^ coskia + {^Ai — — ^ cosA^a = (6.8) 

always holds. In case cos k\ct ^ and cos k 2 a ^ ; then A± = 1 if and only 
if A 2 = 1, and this can only happen ift = and Hex = 0. 
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Proof. First, multiply the equation F^(x,t) = by 1 — ijAie~ lkia , express 
this using the representation (16. 7p and take the imaginary part of the result. 
This leads to the formula 

(Ax + ——] cos k 2 a = — sin(A;2 + ki)a — 7 sin(/c 2 — ki)a. (6.9) 
V A X J 7 

In the same way, multiplying by 1 — i^A 2 e~ %k,ia and subsequently extracting 
the imaginary part yields 

( A 2 + —r-\ cos kia = — sin(/c 2 + ki)a + 7 sin(/c 2 — ki)a. (6.10) 
V AoJ 7 



If cos k 2 a = 0, so that in particular sin k 2 a = ±1, it follows from (16.91) . using 
the formulas for the sine of the sum and difference, that cos k\a = 0. In the 
same way, using (16.101) . if cos k\a = 0, then cos k 2 a = 0. This proves the first 
assertion in the proposition. 

For future reference, notice that it follows by subtracting ( I6.9P from ( I6.10p 
that 

7 sin(A;2 - h)a = - (^A 2 + j cos k x a - - (^A 1 + -|- j cos k 2 a. (6.11) 

Taking the real and imaginary parts of (16. 7p and setting them equal to 
zero yields 

1 — 7A1 sin k±a — r yA 2 sin k 2 a — A\A 2 cos(k\ + k 2 )a = 0, (6-12) 
jAi cos kia + •yA 2 cos k 2 a — A\A 2 sin(fci + k 2 )a = 0. (6.13) 

Multiply the first equation above by sin(fci + k 2 )a, the second equation by 
cos(/ci + k 2 )a and form the difference of the outcomes. The formula 

sin(fci + k 2 )a = •yA 2 cos k\a + 7A1 cos k 2 a (6.14) 

emerges from these machinations. But, (I6.13P implies 

1 1 
sin(A; 1 + k 2 )a = — cos k\a + — cos k 2 a. (6.15) 

A 2 A\ 



The last two equations taken together imply ( 16. 8 p 

This proves the second assertion of the proposition. Finally, it is clear 
that A\ = A 2 — 1 if and only if t = and Rex = 0. □ 
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Proposition 6.4. Suppose that coskia = cosk 2 a = 0. It follows that k\ 
and k 2 are commensurable, that p\ and p 2 are both odd and that a is an odd 
multiple ofnX/2. 

Proof. If cos k\Ct = cos k 2 a = 0, then there exist integers m and n such that 

71 

kia = (2m + 1)"^") 

k 2 a = (2n+l)|. 

Thus it is clear that ki and k 2 are commensurable and 

Pi ki 2m + 1 
p 2 k 2 2n + 1 ' 

whence 

p 1 (2n + 1) =p 2 (2m+ 1). (6.16) 

It follows that p 2 — pi is even, and since they are also relatively prime, they 
both must be odd. Next, equation ( I6.16P implies that p\ = (2m + l)/c, where 

c = gcd(2m + l,2n+ 1). 

Consequently, 

2m + 1 7r 2m + 1 tt\ tt\ 
01 ~ k~i 2 ~ ~p~ x 2 ~ C T' 
which concludes the proof since c is necessarily odd. □ 

Proposition 6.5. Let x(t) be a smooth curve such that F^(x(t),t) = and 
d x F^~(x(t),t) 7^ 0. It follows that lmx'(t) has the same sign as 

(^Ai — — j cos k 2 a = — (^A 2 — — j cos/cia. 

Proof. Let x(t) be a smooth curve such that F^~(x(t),t) = and d x F^~(x(t),t) ^ 
0. Implicit differentiation provides the relation 

d t F+{x(t),t) 



At) 



d x F+(x(t),t) 

kli-fA ie ikia + k%ijA 2 e ik2a - (kf + k^)A 1 A 2 e i( - kl+k ^ a 
k^Aie^ + k 2 i^A 2 e ik * a - {k x + k 2 ) A x A 2 e i( - k ^ k ^ a 

k\~{A x e- ik * a + k^A 2 e~ ikia + i(kf + k\)A x A 2 
k^A^-^ + k 2 <yA 2 e~ ikia + i(h + k 2 )A x A 2 
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&i7v4ie- ife2a + k 2 ^A 2 e- ikia + i(k x + k 2 )A 1 A 2 \ 2 



-ik\a 



+i{k\ + kDA^) (k ll A 1 e tk2a + k 2l A 2 e ik ' a - i(h + k 2 )A 1 A 2 ) 

The imaginary part of the numerator in this last expression is equal to 

kik 2 (k%— kl)^ 2 A±A 2 sm(k 2 —ki)a+kik 2 (ki+k 2 ) 2 AiA 2 (Ai cos k 2 a—A 2 cos kia), 

which is a positive mulitple of, and so has the same sign as, 

7sin(A;2 — ki)a + Ai cos k 2 a — A 2 cos k x a. (6-17) 

Finally, formulas 06. 11 j) and (I6.8P reveal that the expression in (I6.17P is equal 
to 

Ax cos k 2 a — cos k 2 a. 

Ai 

□ 

Corollary 6.1. Suppose either that k\ and k 2 are not commensurable, or 
that they are commensurable with p\ and p 2 having opposite parity. Let x(t) 
be a smooth curve such that F±(x(t),t) = 0. // either t ^ or if Hex ^ 0, 
then lmx'(t) ^ 0. 

Proof. We know from Lemma 13.11 that the roots of F + (-,t) are all simple 
unless k\ and k 2 are commensurable and p\ and p 2 are both odd. It fol- 
lows that the same is true for Thus, it must be the case that 
d x F^~(x(t), t) ytz o. The last three propositions can now be brought to bear 
to establish the claim. □ 

To recapitulate, it has been shown, except in the commensurable case with 
Pi and p 2 both odd, that if either t ^ or if Rex ^ 0, then Imx'(t) ^ 0, so 
the poles of the solution u + are always moving vertically. The same is true in 
the case that p\ and p 2 are both odd, so long as the imaginary part of the pole 
is not an odd multiple of 7rA/2. The same conclusions are true about the poles 
of u~ . The calculations leading to this conclusion, while not the same as those 
for u + , are analogous enough that we pass over the details. Interestingly, in 
the commensurable case when p\ and p 2 have opposite parity, it turns out 
that the poles of u~ and those of u + bear a very simple relationship to one 
another. 
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Proposition 6.6. Suppose that pi and P2 have opposite parity (one odd, one 
even). It follows that there exists fleR such that 

F + (x -iO,t) = F~(x,t) (6.18) 

for all x G C and tel. In other words, the poles of u~ are precisely given 
by a vertical translation of the poles of u + . 

Proof. Suppose 9 G R is such that 

f 1 (x-i9,t) = f 1 (x,t), (6.19) 
h(x-i6,t) = -f 2 {x,t), (6.20) 

for all x G C and tel. It would then follow that (16 . 18[) holds for all x G C 
and tel. The same would be true if we had instead 

f 1 (x-i9,t) = -f 1 (x,t), (6.21) 
f 2 (x-i9,t) = f 2 (x,t). (6.22) 

For (I6.19P and (I6.20p to be valid, it is necessary and sufficient that 

exp(ifci0) = exp(zpi0/A) = 1, (6.23) 
exp(ik 2 9) = exp(ip 2 9/\) = -1. (6.24) 

For these latter conditions to hold, it is necessary for there to be two integers 
m and n such that 

P10/A = 2rmr, (6.25) 
p 2 9/X = (2n + l)vr, (6.26) 

or, what is the same, 

J^ = 2m = 2n + l 
Att p 1 p 2 

If pi is even and p 2 is odd, it is clear that one may choose appropriate values 
of 9, m and n so that the last equation holds 

In the opposite case, if p\ is odd and p 2 is even, a similar argument shows 
that there exists 9, mand n such that (16.211) and (16.221) hold. □ 



28 



It remains to consider the case where k\ and k 2 are commensurable, with 
Pi and P2 odd. To establish the existence of poles with non-trivial vertical 
movement, it suffices by Proposition 16.41 to establish the existence of poles 
whose imaginary parts are not an odd multiple of n\/2 (with either t 7^ or 
with non-zero real part). 

As shown in [3], for the two-soliton solution of the KdV-equation in the 
commensurable case, there are always poles with vertical movement, and at 
least one pole moving precisely horizontally. The proof of this fact requires 
the full force of the delicate and technical analysis in [9J. We would like to 
avoid such technical calculations in this paper to the extent possible. 

Thus, for the modified KdV-equation (I2.ip . we should not expect that 
all the poles will be moving vertically in the remaining cases. The fact that 
in the case of opposite parity, all the poles move vertically (except if t — 
and Rex = 0) is already an interesting difference in behavior between the 
two equations. 

It turns out that in the commensurable case, with p\ and p 2 both odd, 
the movement of the poles of u + can in fact be reduced to the movement of 
the poles of the 2-soliton solution of the KdV-equation if p 2 — P\ G 4N. The 
same is true for u~ if p 2 + P\ G 4N. Here are the precise statements. 

Proposition 6.7. Suppose that pi and p 2 are both odd. If p 2 — p\ G 4N then 
there exist 6\,6 2 G M such that 

F+(x-i9 u t) = l + >yf 1 (x,t) + >yf 2 (x,t)+f 1 (x,t)f 2 (x,t) (6.28) 

and 

F+(x-i9 2 ,t) = l + >yf 1 (x,t) + >yf 2 (x,t)+f 1 (x,t)f 2 (x,t) (6.29) 
for all x G C and tel. 
Proof. Suppose 6\ G 1R is such that 

/i(z-zM) = -z/iOM), (6-30) 
f 2 {x-i6 ll t) = -if 2 {x ) t) ) (6.31) 

for all x G C and tel. It would follow that 06.281) holds for all x G C and 
tel. In addition, if we have 

fi{x-ie 2 ,t)=if 1 (x,t\ (6.32) 
f 2 (x-id 2 ,t)=if 2 (x,t), (6.33) 
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then ( I6.29P would be true. 

For the system (I6.30p - (l6.3ip to be valid, it is necessary and sufficient that 

exp(iki9) = exp(ipi9 / A) = —i, (6.34) 

exp(ik 2 9) = exp(ip 2 9 / A) = —i. (6.35) 

For this, we need to find two integers m and n such that 

^ = (4m - (6.36) 

P 4 = (4n - D|, (6.37) 

which is the same as asking for two integers m and n such that 

29 Am -I An -I 



(6.38) 



Att pi p 2 
For the latter to hold true, it must be the case that 

4(p2m - pi " } = 1. (6.39) 

P2 ~Pl 

Since p\ and p 2 are relatively prime, there exist integers r and s such that 

rp 2 + spi = 1; 

hence, simply take 

G°2 - Pi)r 



m 



n 



A 

(p 2 -pi)s 
A 



The proof for F 2 + is similar, but with Am + 1 replacing Am — 1 and An + 1 
replacing An — 1. □ 

Proposition 6.8. Suppose that pi and p 2 are both odd. Ifp 2 + Pi € 4N i/ien 
i/iere exisi 9\,9 2 £R swc/i i/iai 

Ff(x - i9 1 , t) = 1 + 7/i(z, *) + 7/ 2 (x, t) + t)/ 2 (x, t) (6.40) 

and 

F~(x - i0 2 , t) = 1 + 7/1 (z, *) + 7/2(2:, t) + /1 (x, t)/ 2 (x, t) (6.41) 
for all x G C and tel. 
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The expression on the right side of the four formulas (I6.28P , ( I6.29P , (I6.40p , 
and (KT4"T|) . i.e. 



I+7/1OM) + jf 2 (x,t) + h{x,t)f 2 (x,t), (6.42) 

is exactly the function F in formula (2.13) of [9] whose zeros correspond to 
the poles of the 2-soliton solution of the KdV-equation. Thus, we may use 
the results of [9] to describe the behavior of the poles of in the cases under 
consideration. 

More precisely, we may now affirm that in the case where p\ and p 2 are 
both odd and p 2 —pi G 4N, the solution u + has two poles moving horizontally 
on the line Imx = ttX/2 and also on the line Imx = —nX/2. There are 
2(pi + p 2 — 2) other poles with imaginary part between ±7rA, and they will 
move vertically for all t ^ 0. The same will be true for u~ in the case p\ and 
p 2 are both odd if p 2 + P\ G 4N. These configurations will repeat with 2irXi 
periodicity, so that horizontally moving poles are found with imaginary part 
equal to every odd multiple of nX/2. 

The last case is one where four poles meet at t = 0, as described in the 
previous section, at poles on the imaginary axis whose imaginary part is an 
odd multiple of ttX/2. At each of these points, the analysis shows that (at 
least) two of these poles have vertical movement. We refrain from going 
into the detailed considerations needed to establish that all the poles move 
vertically for all t ^ 0, excepting the two poles approaching the singular 
points horizontally. In any event, we have shown in this case the existence 
of poles with vertical movement. 

7 Finite time blowup of solutions 

The results of the previous section immediately imply that there exist complex- 
valued solutions to (12.1 p on R which blow up in finite time. Indeed, fix a G R, 
and let u be given by 

u(x, t) = m ± (x — ia, t), (7.1) 

where we may use either of the two soliton solutions u + or u~ . As long as 
the set {(x — ia,to) : x G M} does not contain a pole of m ± , u is a smooth, 
complex- valued solution of (12. ip for t in a neighborhood of to- Since at any 
given time, the imaginary parts of the collection of all poles of u + or u~ form 
a discrete set, and since both u + and u~ have poles whose imaginary parts 
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move continuously in time, there exist a and to such that [x — ia, to) is a 
pole of u + , say, but (x — ia,t) is not a pole of u + if t is sufficiently close 
to, but not equal to, to- If follows that the solution u defined in (17. ip with 
this choice of a is a regular complex-valued solution of (12.11) on K, decaying 
exponentially to zero as x — > ±00 for t close to to, but which is singular at 
t = t . In other words, the solution blows up in finite time. 

It is interesting to note that this result of singularity formation for complex- 
valued solutions of the mKdV-equation can be interpreted as a blow-up re- 
sult for real-valued solutions of a system of dispersive equations. Let u be 
a solution of (12. ip and let r = Reu and s = Imu. It follows that r and s 
satisfy the real-valued system 

r t + r xxx + 6(r 2 - s 2 )r x - 2rss x = 0, (7.2) 
St + s xxx + 2rsr x + 6(r 2 - s 2 )s x = 0. (7.3) 

Thus, we have shown that this coupled, dispersive system admits real-valued 
solutions (exponentially decaying in space) which blow up in finite time. 

8 Some formal calculations 

(All the computations in this section have been carried out using MAPLE.) 
As noted at the end of Section [2J the interaction time for the two-soliton 
solutions m ± , given in (I2.13P and (I2.14p . is t = 0, the center of the interaction 
is x = and at t = 0, the solution is even in x. The explicit formulas for 
m ± allow us to observe and calculate certain aspects of these solutions at the 
moment of interaction. In particular, it is interesting to know whether there 
is a single maxima during the interaction or not, and it is likewise interesting 
to know the speed of the two solitons at the moment of interaction. 

In the case of u + , one observes that that the solution has one centered 
maximum if the ratio k 2 /ki is large enough (bigger than around 2.6), and 
two symmetrically located maxima for smaller values of ^2/^1- Partial confir- 
mation of this can be obtained by computing u xx (0, 0). Since, by symmetry, 
«+((), 0) = 0, the sign of the second derivative will tell us if it is a local 
maximum or a local minimum. A MAPLE-implemented calculation shows 
that 

*4(0, 0) = -(k 2 - k^kl - 3hk 2 + k 2 2 ). 

Therefore, if 1 < k 2 /k 1 < (3 + y/5)/2, then u+.(0,0) > 0, which means that 
u + (-, 0) has a local minimum at x — 0. Thus, there are (at least) two maxima 
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at the moment of interaction. If k 2 /k\ > (3 + v5)/2, then 0) has a local 
maximum at x — 0, which is consistent with there being a single maximum 
at the moment of interaction. 
Similarly, we calculate that 



■ft 



(0, 0) = -(k 2 + h){kl + 3hk 2 + k 2 2 ) 



which means it~(-, 0) has a local maximum at x = for all values < k\ < k 2 , 
which is consistent with the graphical observations of the solution itself. 

If we are interested in the "speed" of the two-soliton solution at the 
moment of interaction, one idea is to calculate the speed of the maximum 
or the local minimum of the solution as (x,t) approaches (0,0). One might 
argue that the movement of the maximum is some kind of speed. It must 
be acknowledged that the interpretation is less clear when one is tracking a 
local minimum. To calculate this speed, suppose y(t) = y±(t) is a real-valued 
curve such that u x (y(t),t) = 0, i.e. y±(t) is always at an external point of 
the solution u = u ± . Suppose also that y±(0) = 0, which is to say the curve 
is at the interaction center at the interaction time t = 0. Differentiating with 
respect to t, we see that u xx (y(t),t)y'(t) + u xt (y(t),t) = 0, or 



HQ = u xt (y(t),t) 
' ' u xx (y(t),t)' 



This gives, in turn, 



,'(0) = (8.1) 
u xx {0,0) 

which might be thought of as representing the speed of the two-soliton solu- 
tion at the moment of interaction of the two solitons. The value of y'(0) can 
be calculated explicitly from (I2.13P and ( 12.14]) . The results are as follows, 
for both it + and u~, 



<4(o) 
y-(o) 



k\ - 3hk 2 + k\ 
k\ + 3hk 2 + k\ 



For it - , where there is always a maximum at x = at the moment of in- 
teraction, the maximum is moving with a posistive speed. What that speed 
represents is not entirely clear. In the case of u + , where the midpoint is a 
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local maximum only if k 2 /k 1 > (3 + Vb)/2, we see that for these values of 
< ki < k 2 , we have indeed y' + (0) > 0. On the other hand, y' + (0) < for at 
least some values of < ki < k 2 with k 2 /ki < (3 + \/5)/2. (The lower bound 
on k 2 /ki for which this speed is negative is around 2.15.) This negative speed 
represents the speed of the local minimum, between the two maxima. It is 
curious that in some cases this minimum is moving backwards. 

It is also interesting to do this for the two-soliton solution of the KdV- 
equation, a calculation which was not carried out in [9]. In this case, we 
obtain 

. fig i- ^i^ 2 fc 2 

y ~~ op _ p 

As is known, the two-soliton solution of the KdV-equation has one maximum 
at the interaction time if k 2 /ki > V3 and two maxima, symmetrically located 
about the interaction center, if 1 < k 2 /k\ < y/3. Here we see that y'(0) < 
for a/ 1 + y/2 < k 2 /k\ < y/3. In these case, of course, it is the minimum 
which is moving backward. 
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